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I.  INTRODUCTION 

1.1.  Statement  of  the  Problem 

The  procedure  of  blood  flow  measurement  is  well  established  and  extends  into  many  areas 
of  medical  practice.  Direct  indication  of  the  patient's  hemodynamic  status  (and,  less  directly,  the 
overall  status)  results  from  analysis  of  basic  blood  flow  characteristics:  velocity,  pulsatility, 
resistance  and  so  on.  Vascular  disorders  such  as  plaque  accumulation,  stenosis  or  aneurysm  can 
be  diagnosed  and  monitored  by  imaging  or  otherwise  evaluating  blood  flow  characteristics. 
Quantification  of  large  vessel  blood  flow  is  integral  in  determining  the  course  of  treatment  for  a 
specific  disorder  or  patient.  Current  techniques  used  to  obtain  quantified  flow  data,  however,  are 
limited  to  larger  blood  vessels  systems  such  as  the  carotids,  femorals,  aorta,  renal  and  hepatic 
vessels,  and  are  often  restricted  to  the  more  superficially  located  vessels  found  easily  by  clinical 
exam  (ie,  feeling  for  the  pulse). 

Blood  perfusion  of  tissue,  that  is,  blood  flow  through  the  tiny  microcirculation,  is  an 
important  area  of  clinical  investigation.  The  viability  of  a  region  of  tissue  is  directly  related  to  the 
perfusion  of  (and,  thus,  the  oxygen  supply  available  to)  that  region.  For  example,  trauma  to  one 
of  the  large  distributing  vessels  to  a  distal  limb  can  result  in  underperfusion  of  the  tissue  in  that 
region,  resulting  in  ischemia  and  eventual  necrosis  necessitating  amputation.  Similarly,  chronic 
occlusive  coronary  disease  results  in  underperfusion  of  a  region  of  myocardium.  Early  detection 
of  the  low  perfusion  could  avoid  life-endangering  complications  of  myocardial  infarction  through 
appropriate  therapy.  Another  example  is  the  neovascularity  characteristic  of  many  tumors.  This 
new  vasculature  supplying  nutrients  to  the  rapidly  growing  tissue  comprises  mostly  small  vessels 
which  elude  detection  by  conventional  means.  Unfortunately,  available  methods  for  evaluating 
blood  perfusion  lack  the  highly  desirable  characteristics  of  noninvasiveness  and  quantifiability. 

An  ideal  blood  perfusion  diagnostic  device  for  use  in  triage  situations  would  be 
uncomplicated  to  use  under  high  stress  and  often  mobile  conditions.  Such  a  device  would 
noninvasively  supply  information  about  the  hemodynamic  status  of  deep  tissue.  Accurate  and 
precise  measures  would  be  obtained  which  would  be  immediately  recognizable  by,  and  significant 
to,  trained  personnel  as  reliable  indicators  of  tissue  survivability.  The  device  would  also  be 
versatile  enough  to  remain  with  a  patient  for  extended  monitoring  and  should  be  compatible  with 
commonly  used  display  and  data  storage  devices.  Currently  no  system  or  device  meets  these 
requirements. 

Until  roughly  the  last  decade,  limitations  in  blood  flow  diagnostic  technologies  have  made 
analysis  of  blood  perfusion  (flow  in  a  bed  of  the  smallest  vessels:  arterioles,  capillaries,  venules) 
clinically  impractical.  Other  noninvasive  techniques  addressed  the  desire  to  obtain  absolute  flow 
information  from  the  microcirculation  but  were  limited  in  robustness  (e.g.,  plethysmography, 
which  required  a  stable  hemodynamic  state,  an  appreciable  amount  of  time  to  generate  the  data,  and 
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awkward,  bulky,  somewhat  fragile  equipment),  reliability  (e.g.,  133Xenon  clearance)  and  regions 
of  measurement  (e.g.,  thermal  clearance  and  laser  Doppler  flowmetry  which  measure  superficial 
flow). 

Limited  progress  in  noninvasively  measuring  microcirculatory  flow  via  ultrasound  has  been 
reported.  The  obstacles  in  the  way  of  such  a  diagnostic  procedure  derive  from  the  physiology  of 
the  measured  system:  slow  blood  velocity  and  low  volume.  Capillary  blood  velocity  is, 
characteristically,  no  greater  than  about  2  mm/sec;  the  relative  blood  volume  of  tissue  at  resting 
perfusion  is  about  l/1000th  that  of  tissue.  A  Doppler  system  operating  at  5  MHz  would  produce 
frequency  deviations  on  the  order  of  tens  of  Hertz  for  velocities  on  the  order  of  1  mm/sec,  below 
the  resolution  of  conventional  systems.  Also,  tissue  movement  from  either  gross  patient  motion  or 
from  autonomic  vasomotive  processes  can  "bury"  the  signal  of  interest  in  artifact,  particularly 
because  of  the  small  capillary  blood  volume  in  the  measurement  region  (which  produces  a  much 
decreased  blood-to-tissue  echo  ratio). 

Some  preliminary  theoretical  work  (Bums  and  Reid,  1989)  has  addressed  the  small 
frequency  shift  problem  via  stable  oscillators.  The  proposed  system  was  expected  to  resolve  shifts 
of  less  than  10  Hertz  from  a  1  MHz  carrier;  development  and  tests  of  the  system  have  not  been 
published.  Dymling,  et  al  (1991),  used  a  continuous  wave  Doppler  system  without  the  customary 
low-pass  filter  and  reported  promising  results.  The  method  is  limited,  however,  by  tissue  motion 
artifact. 

Use  of  a  contrast  agent  to  enhance  the  signal- to-clutter  ratio  of  capillary  blood  flow 
addresses  the  small  signal  problem  inherent  in  attempting  to  measure  microcirculatory  perfusion  via 
Doppler  ultrasound.  Reports  by  Monaghan,  et  al  (1988),  Ong,  et  al  (1984),  Ten  Cate,  et  al  (1984) 
have  at  least  demonstrated  the  feasibility  of  contrast  enhancement  in  microcirculatory  flow 
measurement.  However,  the  small  blood  volume  involved,  along  with  tissue  motion  artifact,  can 
confound  the  frequency  shift  or  intensity  of  the  Doppler  signal  and  still  lead  to  ambiguous  results. 

Another  limitation  of  existing  ultrasound  techniques  is  lack  of  consistency  in  the 
measurement  system.  Changes  in  the  quality  of  coupling  between  the  transducer  and  skin  or  in  the 
angle  of  the  incident  beam  can  cause  variations  in  the  incident  and  reflected  ultrasonic  energy.  This 
inconsistency  is  inherent  in  all  of  the  present  systems,  including  those  utilizing  contrast 
enhancement.  Therefore,  absolute  quantification  of  signal  intensity  and,  thus,  absolute 
quantification  of  perfusion,  is  not  possible  with  current  methods. 

Analysis  by  Schrope  and  Newhouse  (1993)  demonstrated  that  intervening  tissue  produced 
negligible  nonlinearity  in  the  propagated  signal  as  compared  with  a  selected  contrast  agent  in  the 
capillaries.  The  nonlinear  contrast  agent  generates  a  second  harmonic  component  which  allows 
elimination  of  the  clutter  signal  returned  at  the  fundamental  frequency  by  vasomotive  processes  and 
intervening  tissue.  Utilizing  this  nonlinear  contrast  agent  and  recently  developed  algorithms,  our 
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proposed  system  overcomes  the  obstacles  which  hindered  previous  methods  in  order  to  provide  a 
technique  for  the  unambiguous  quantification  of  tissue  blood  flow  measurement. 

The  research  covered  by  this  report  is  directed  to  developing  a  device  which  will  provide  an 
unambiguous  quantitative  measurement  of  tissue  perfusion  and  meet  the  other  requirements  as 
stated  in  the  preceding  section  via  Doppler  ultrasound  techniques.  Previous  work  (Schrope,  et  al, 
1992;  Schrope  and  Newhouse,  1993)  demonstrated  the  advantage  of  using  a  nonlinear  contrast 
agent  and  analyzing  the  resulting  second  harmonic  component  of  the  return  signal  (echo).  The 
analysis  technique  greatly  reduces  the  "clutter"  produced  by  echoes  from  surrounding  tissue  and 
vasomotive  processes  providing  a  well-differentiated  signal  due  only  to  blood  flow.  Additionally, 
a  unique  approach  to  the  problem  of  quantifying  the  signal  intensity  at  the  target  tissue  and,  thus 
the  absolute  measure  of  blood  flow,  will  result  in  a  calibrated  system  which  offers  useful  and 
reliable  quantitative  information  about  deep  tissue  blood  flow. 

1.2.  Phase  I  Goals 

As  delineated  in  the  Phase  I  Proposal,  our  method  of  investigation  in  this  project 
encompassed  four  technical  objectives: 

1.  Assess  the  appropriateness  and  technical  feasibility  of  the  second  harmonic  ultrasonic 
microcirculatory  blood  flow  measurement  in  the  application  of  trauma  medicine. 

To  achieve  this  goal  we  began  with  a  detailed  discussion  of  the  "Soldier  Individual 
Computer"  project  with  a  member  of  the  team  in  the  Army.  We  learned  much  more  about  the 
ultimate  goals  of  this  large  project  than  were  available  in  the  RFP,  and  we  were  able  to  adjust  our 
priorities  and  designs  accordingly. 

2.  Refine  and  adapt  algorithms  for  second  harmonic  tissue  perfusion  measurement  to  be  compatible 
with  commonly  available  display  and  data  storage  devices,  such  as  used  by  the  military  and 
commercial  sectors. 

Again,  the  discussion  with  the  Army  project  representative  clarified  this  goal  for  the 
project.  Specifically,  we  learned  that  the  Army  desired  a  specific  data  output  for  integration  into 
the  Soldier  Individual  Computer  system.  We  endeavored  to  ultimately  provide  this  output  from 
our  algorithms. 

3.  Build  prototype  device  that  is  economical,  durable,  easy-to-use  and  compatible  with  existing 
systems. 

As  stated  in  the  proposal,  our  phase  I  effort  was  not  directed  at  "reinventing  the  wheel"  - 
hence,  we  chose  mostly  off-the-shelf  components  for  our  feasibility  testing  system.  Here  we 
report  several  designs  for  systems  that  were  conceived  during  the  project  period;  in  our 
experimental  methods  section  we  describe  the  design  actually  used  for  the  experimentation  and 
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testing.  Eventually  our  technology  and  algorithms  will  be  integrated  into  a  small  DSP  computer 
chip  and  an  external  sensor,  for  integration  into  the  Soldier  Individual  Computer. 

4.  Test  prototype  on  in  vitro  perfusion  models. 

This  project  entailed  extensive  testing  of  the  concept  under  question,  that  is,  the  second 
harmonic  perfusion  measurement.  As  will  be  described  later  in  this  report,  we  endeavored  to 
carefully  determine  the  useful  limits  of  the  system,  in  terms  of  frequency,  contrast  agent 
concentration,  input  power,  etc.. 

H.  SECOND  HARMONIC  SYSTEM  DESIGN 

In  this  section  is  a  summary  of  several  designs  for  a  system  to  measure  and  quantify  blood 
perfusion  using  the  second  harmonic  signal  detected  from  the  region  of  interest.  In  general,  we 
proposed  three  different  types  of  systems. 

Our  first  system,  schematically  demonstrated  in  Figure  1,  is  based  on  apparatus  currently 
used  in  the  laboratory  with  the  addition  of  a  few  enhancements  such  as  PC  interface,  improved 
precision  equipment,  etc..  The  advantages  of  such  a  system  are  that:  it  requires  minimal 
engineering  effort,  leaving  more  time  available  for  extensive  experimentation  and  understanding  of 
the  contrast  agent;  it  is  very  flexible  due  to  the  general  nature  of  the  lab  test  equipment;  it's  A/D 
conversion  and  PC  interface  allows  storage  of  raw  data  which  can  be  saved  and  manipulated  in  the 
future  as  new  and  improved  algorithms  become  available.  The  main  disadvantage  of  this  system, 
also  related  to  its  flexibility  and  general  nature  as  test  equipment,  is  its  high  cost. 

The  system  depicted  in  Figure  2  is  actually  a  power  estimator  of  both  the  first  and  second 
harmonic.  Again  we  convert  our  data  to  digital  format  for  analysis;  however,  the  A/D  converters 
are  employed  at  a  much  reduced  speed  of  less  than  1  MHz  and  8  bit  A/D  conversion  (versus  at  least 
12  bit  needed  for  the  first  approach).  The  main  advantages  of  this  system  are:  it  is  truly  a  prototype 
system  that  can  lead  to  product  design  if  the  approach  is  found  successful;  it  retains  the  feature  of 
the  first  approach  in  that  the  first  and  second  harmonic  amplitudes  can  be  estimated  from  within  a 
specified  range  cell;  integration  of  other  data  analysis  features  (for  example,  Doppler  or  the 
attenuation  analysis  to  be  discussed  later  in  this  report)  is  rather  straightforward.  However,  this 
system  requires  significantly  more  engineering  effort,  which  may  not  be  appropriate  at  this  stage  of 
the  research.  Also,  this  design  is  less  flexible  than  the  first  (mainly  in  controlling  the 
characteristics  of  the  transmitted  signal  such  as  frequency,  burst  length,  power,  etc.).  Of  course, 
these  parameters  could  be  altered  by  modifying  the  system,  but  in  the  first  design  the  parameters 
could  be  altered  continuously  and  within  the  same  experiment. 

Referring  again  to  Figure  2,  the  details  of  the  design  are  as  follows.  The  transmitter 
transmits  bursts  at  co  =  2icf  at  a  given  repetition  rate  somewhere  between  500  and  2000  Hz.  The 
received  signal  is  amplified  (amplifiers  not  shown)  and  split  into  two  signals.  In  the  first  branch 
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(upper)  the  fundamental  harmonic  signal  is  processed,  and  in  the  second  branch  (lower)  the  second 
harmonic  signal  is  processed.  The  signals  are  processed  in  quadrature,  so  restoration  of  the 
signals’  envelope  and  with  it  the  power  are  possible.  We  in  fact  created  a  power  meter  for  the 
fundamental  and  second  harmonic  signals.  In  order  to  restore  the  signal  power,  we  take  the 
signals  from  the  quadrature  outputs  and,  using  a  lookup  table,  estimate  the  returned  signal 
envelopes  of  the  gate  signals.  Averaging  of  consecutive  returned  signals  improves  the  robustness 
of  the  results. 

Theoretically,  the  development  is  as  follows.  The  returned  signal  within  the  gate  region 
can  be  described  as: 

r(t)  =  A(t)  sin  (cot  +  0)  (II.  1) 

After  multiplication  with  the  two  quadrature  signals  sin  (cot),  cos(cot),  we  obtain  the  following 
signals  (after  the  mixer  and  low  pass  filters): 

r  i(t)  =  A(t)sin((J>)  (II.2) 

r2(t)  =  A(t)cos((J))  (II. 3) 

Restoration  of  A(t)  is  easily  achievable  through 

A(t)  =  (r12(t)  +  r22(t))0'5  (11*4) 

which  can  best  be  accomplished  via  lookup  tables.  From  A(t)  we  estimate  its  average  within  the 
gate.  The  same  analysis  applies  to  the  second  harmonic  signals.  Note  that  the  0.5  MHz  filters  are 
anti-aliasing  filters  and  their  function  is  to  retain  the  original  information  contained  in  the  returned 
signals  for  the  fundamental  harmonic  in  the  upper  branch  and  the  second  harmonic  in  the  lower 
branch.  The  bandwidth  of  the  signals  is  mainly  defined  by  the  burst  length.  For  example,  for  a  3 
MHz  signal  with  12  cycles  burst  length,  the  bandwidth  after  down-conversion  will  be  for 
approximately  3/24  MHz,  which  is  appropriately  less  than  the  250  kHz  low  pass  filter  suggested  in 
the  design.  Even  for  cases  when  the  returned  signal  would  require  larger  bandwidths,  the  LPF 
would  just  cut  the  average  power  in  a  predictable  manner  and  accurate  estimates  of  the  power  of  the 
fundamental  and  second  harmonic  signals  would  be  possible.  (Note  that  we  are  interested  not  in 
the  spectrum  shape  but  in  its  power.)  The  above  design  assumes  that  we  will  need  both 
fundamental  and  second  harmonic  signals  for  quantification,  since  at  this  time  our  theory  utilizes 
the  ratio  of  the  two. 

As  mentioned  above,  we  have  reduced  significantly  the  requirements  (and  thus  the  cost  and 
fragility)  for  the  A/D  system  -  from  sampling  rates  of  up  to  50  MHz  and  12  bit  A/D  conversion  we 
can  come  down  to  rates  of  1  MHz  and  8  bit  A/D.  We  also  eliminated  the  FFT  and  thus  eliminated 
its  computational  time,  space,  and  uncertainty.  However,  we  will  need  a  total  of  four  A/D 
converters  (but  since  the  cost  of  a  smaller  bit  converter  is  less,  this  should  be  more  than 
compensated  for).  We  could  even  reduce  this  number  to  two  by  switching  a  pair  of  them  from  the 
fundamental  harmonic  to  the  second  harmonic  every  burst,  so  that  half  of  the  returned  signals  are 
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processed  on  the  fundamental  harmonic  branch  and  half  on  the  second.  This  will  double  the  data 
acquisition  time,  however. 

The  final  system  is  shown  in  Figure  3.  This  one  is  a  Doppler  system  that  is  capable  of 
providing  information  regarding  the  contrast  agent  flow  in  addition  to  its  fundamental  and  second 
harmonic  backscatter  and  attenuation  characteristics.  Its  implementation  is  simpler  in  that  the  A/D 
conversion  is  done  at  very  low  speed  (on  the  order  of  kHz,  compared  to  1  MHz  for  the  second 
approach  and  50  MHz  for  the  first  approach)  and  only  8  bits  are  necessary.  Also,  if  desired,  it 
would  be  simple  to  piggy-back  it  onto  the  second  system  if  desired.  However,  Doppler  detection 
of  low  velocity  flows  is  still  not  definitively  proven  (although  use  of  the  second  harmonic  should 
circumvent  several  of  the  problems  associated  with  it). 

In  Figure  2,  the  points  labelled  (A)  and  (B)  are  shown  as  start  points  for  the  design  in 
Figure  3.  All  of  the  fundamental  harmonic  received  signal  processing  (the  upper  two  mixers)  can 
be  eliminated  if  we  design  the  system  only  for  Doppler.  Also  the  high  speed  A/D  can  be  removed 
from  the  design,  replaced  by  much  lower  speed  A/D  converters.  Note  that  this  design  is  also  a 
quadrature  design,  and,  as  a  Doppler  system,  provide  information  regarding  approaching  and 
receding  flow.  If  no  directionality  is  needed,  the  design  becomes  that  much  simpler;  only  branch  A 
and  no  90  degree  phase  shift  is  needed.  Also,  only  one  A/D  is  necessary. 

m.  EXPERIMENTAL  METHODS 
DLL  Equipment 

The  components  which  constituted  the  system  varied  according  to  the  system  configuration 
under  test  for  an  experiment.  The  description  of  the  components  and  the  function(s)  provided  by 
each  follows: 

LeCroy  9410  oscilloscope:  4-trace,  150  MHz  DSO.  The  waveform  function  option  provided  real 
time  FFT  (amplitude  and  power)  as  well  as  averaging,  among  others.  User-controlled  cursors 
allowed  accurate  and  repeatable  measurements  directly  from  the  'scope  display.  An  external 
trigger  input  provided  capability  to  synchronize  the  sweep  to  "difficult"  signals. 
Hewlett-Packard  81 16A  function  generator:  50  MHz,  programmable  function.  The  81 16A 
provided  both  continuous-wave  (CW)  and  tone  burst  capability.  Frequency,  amplitude, 
number  of  cycles,  burst  repetition  rate  and  other  parameters  were  user-controlled.  A  gate 
output  provided  a  TTL-level  signal  synchronized  with  the  tone  burst:  the  gate  rising  edge 
occurred  at  the  start  of  the  burst;  the  falling  edge  occurred  at  time  =  one-half  of  the  burst 
repetition  interval. 

DC  power  supply:  a  standard  bench  supply;  triple-output  (0  to  6,  0  to  -25,  0  to  25). 

Preamplifiers:  mini-circuits  ZFL-500LN  and  Perry  Amplifier  070/40.  The  preamplifiers  were  used 
either  singly  or  cascaded  to  provide  19,  40,  or  70  dB  gain. 
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Transmit/receive  switch:  fabricated  in  house.  The  switch  was  used  in  single  transducer 

experiments  to  isolate  the  receiver  (either  a  preamplifier  or  the  oscilloscope)  from  the  transmit 
burst.  Isolation  was  approximately  60  dB;  the  gate  output  from  the  function  generator 
provided  the  control  signal  to  the  switch.  The  duration  in  which  the  transmit  signal  could  pass 
through  the  switch  to  the  transducer  was  adjustable  from  approximately  2  pis  to  1  ms. 

Test  tank:  fabricated  in  house.  This  component  provided  a  test  environment  for  the  bench-top 
experiments  involving  equipment  verification,  contrast  agent  characterization  and  preliminary 
feasibility  information  regarding  the  various  system  configurations  and  contrast  agent 
preparations.  The  tank  held  two  transducers  at  a  90  degree  included  angle;  contained  volume 
approximately  equaled  70  ml.  Material:  Plexiglas  and  PVC  fittings. 

Transducers:  Panametrics  V323,  2.25  MHz;  V384,  3.5  MHz;  V310,  5.0  MHz;  V312,  10.0  MHz 
0.25-in.  element,  unfocused  immersion;  V313,  15  MHz  0.25-in.  element,  spherical  focus 
(0.75-in.  focus)  immersion  transducers. 

Contrast  agent:  Schering  AG,  SH  U  508a  (also  known  by  the  trade  name  Levovist),  a  preparation 
consisting  of  granules  of  galactose  and  trace  palmitic  acid,  which  is  prepared  for  its  labelled  use 
by  mixing  a  given  quantity  of  the  dry  portion  with  a  solvent,  either  water,  saline,  or  buffer 
solution. 

Tyrode  buffer  solution:  buffer  used  in  contrast  agent  preparations  for  in  vitro  experiments  (8  g 
NaCl,  1  g  NaHC03,  0.0576  g  NaH2P04  per  L  deionized  H20). 

As  indicated  above,  discussion  with  the  technical  administrator  for  the  grant  offered  further 
understanding  of  the  relation  of  this  project  to  the  umbrella  project  of  the  Soldier  Individual 
Computer  system  with  medical  diagnostic  capability.  This  brought  with  it  the  realization  that  we 
need  not  consider  display  and  storage  systems  and  need  only  supply  a  front-end  system:  signal 
conditioning  for  the  sensor(s)  and  signal  processing  which  would  provide  a  calibrated  output 
directly  related  to  blood  flow  in  the  region  of  tissue  under  measurement.  The  components  listed 
above  (along  with  a  desktop  PC)  served  these  functions  in  this  proof-of-concept  stage. 

III.2.  Description  of  Experiments 

In  general,  the  goal  of  the  experiments  was  twofold.  One,  to  explore  the  limits  of 
usefulness  of  the  experimental  parameters  of  contrast  agent  concentration,  drive  voltage,  incident 
frequency,  and  tone  burst  length  (where  appropriate).  This  was  determined  by  measuring  the 
fundamental  and  second  harmonic  responses  as  well  as  the  mathematically  derived  ratio  of  the  two. 
Once  these  limits  were  determined,  then  we  proceeded  to  the  second  goal,  which  was  to  verify  the 
results  in  a  realistic  ex  vivo  preparation,  where  real  animal  tissues  added  another  level  of 
complexity  to  the  apparatus,  contributing  further  to  attenuation  and  nonlinear  signal  generation. 
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III.2.1.  Contrast  agent  lifetime 

In  order  to  have  confidence  that  recorded  data  reflected  activity  of  controlled  parameters  and 
not  from  unknown  behavior  of  the  contrast  agent  preparation,  experiments  which  explored  the  time 
dependent  characteristics  of  the  contrast  agent  in  solution  were  performed.  These  experiments 
provided  information  which  allowed  us  to  determine  how  long  an  individual  test  trial  could  be,  in 
subsequent  experiments,  before  the  contrast  agent  degradation  significantly  affected  the  results. 

The  general  equipment  configuration  is  shown  in  Figures  4  and  5.  Figure  4  depicts  the 
connections  used  in  a  two  transducer  tone  burst  setup;  Figure  5  shows  the  setup  used  for  one 
transducer.  A  continuous  wave  configuration  was  also  used  and  differed  only  slightly  from  that 
shown  in  Figure  4  in  that  the  transmit  transducer  was  driven  with  a  continuous  sinusoidal  voltage. 

Two  transducers  were  installed  in  the  test  tank.  Figure  6  illustrates  the  holding  fixture  for  a 
transducer.  The  transducers  were  friction-fitted  with  a  rubber  cone  washer  which  was  sized  for  a 
friction  fit  in  the  bore  of  the  holding  fixture.  The  insertion  collar  was  sized  to  be  nearly  an 
interference  fit  with  both  the  bore  of  the  holding  fixture  and  the  body  of  the  transducer,  but  to  still 
allow  easy  insertion  and  withdrawal.  Use  of  the  insertion  collar  ensured  axial  alignment  between 
the  transducer  and  the  holding  fixture  and  the  fixture  arrangement  allowed  a  repeatable  fixation  of 
path  length. 

Prior  to  transducer  insertion,  the  bore  of  the  holding  fixture  was  filled  with  a  coupling 
medium  (degassed,  distilled  water).  Coupling  fluid  displaced  during  transducer  insertion  simply 
drained  from  the  relief  port  in  the  holding  fixture.  Provision  for  a  relief  port  was  necessary 
because  of  the  plastic  film  membrane  which  was  used  early  in  the  project  to  separate  the  coupling 
medium  from  the  solution  under  experiment  in  the  test  tank.  A  piece  of  the  film  was  affixed  to  the 
end  of  the  holding  fixture  with  rubber  cement;  after  sufficient  curing  time,  mild  heat  applied  to  the 
film  shrank  it  enough  to  create  a  taut,  wrinkle-free  septum.  In  later  experiments  we  eliminated  the 
coupling  medium  and  the  membrane  to  allow  direct  contact  of  the  transducer  with  the  test  solution. 

The  transducer  mounting  fixtures  ensured  perpendicular  alignment  of  the  acoustic  axes  to 
avoid  direct  impingement  of  the  fundamental  frequency  energy  from  the  transmit  transducer  onto 
the  receive  transducer.  Acoustic  reflections  from  bubbly  solutions  such  as  the  contrast  agent  occur 
diffusely  (in  all  directions),  so  there  is  no  detriment  to  receiving  echoes  from  a  perpendicular  angle 
to  the  transmitted  energy.  Second  harmonic  signals  would  be  observed  from  a  region  in  the  test 
solution  defined  by  the  intersection  of  the  transducer  beam  patterns.  Various  concentrations  of  the 
contrast  agent  in  Tyrode  buffer  solution  constituted  the  test  solutions  used  in  the  tank. 

Transducer  pairs,  in  the  two  transducer  configuration,  always  differed  in  center  frequency; 
the  lower  frequency  transducer  being  the  transmitter  and  the  higher  frequency  transducer  acted  as 
the  receiver.  The  center  frequencies  and  response  characteristics  of  the  transducers  and  function 
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generator  drive  frequencies  were  chosen  to  minimize  direct  second  harmonic  energy  from  the 
transmitter  and  sensitivity  of  the  receiver  to  the  fundamental  frequency  energy.  This  arrangement 
lent  assurance  that  second  harmonic  energy  was  due  primarily  to  the  contrast  agent  nonlinearity, 
and  that  the  dynamic  range  of  the  acquired  signals  would  not  tax  the  acquisition  equipment. 

The  HP  function  generator  supplied  the  drive  voltage  for  the  transmit  transducer.  In  the 
tone  burst  configuration,  the  function  generator  gate  output  signal  provided  oscilloscope  sweep 
synchronization  to  the  transmit  signal,  which  was  monitored  in  tone  burst  setup  on  the 
oscilloscope.  Continuous  wave  configurations,  of  course,  did  not  require  sweep  synchronization. 

Preamplifiers  provided  the  sensitivity  needed  by  the  system  to  acquire  the  low  level  (on  the 
order  of  tens  of  microvolts,  in  some  experiments)  received  signals.  Use  of  the  two  available 
preamplifiers  allowed  a  selection  of  gains  determined  by  which  preamplifier  was  used  singly,  or  by 
cascading  the  two  preamplifiers.  The  bench  top  DC  power  supply  provided  operating  voltages 
required  by  the  preamplifier.  The  amplified  receive  signal  passed  on  to  the  oscilloscope. 

The  waveform  functions  available  in  the  oscilloscope  allowed  monitoring  of  the  time 
domain  signal  as  well  as  acquisition  of  a  processed  signal  from  which  the  data  were  taken.  The 
built  in  cursor  functions  provided  a  means  of  obtaining  accurate  readings  in  a  repeatable  manner 

An  experimental  trial  consisted  of  acquiring  data  for  a  given  length  of  time  for  a  particular 
set  of  parameters  .  A  test  solution  was  prepared  by  weighing  the  desired  amount  of  dry  agent 
using  a  precision  microbalance  (Sartorius)  and  mixing  this  by  gentle  agitation  with  Tyrode  solution 
(described  above).  The  mixture  was  gently  stirred  until  no  granules  were  visible  and  then  slowly 
poured  from  the  beaker  into  the  test  tank.  To  avoid  reading  the  initial  highly  kinetic  transient 
activity  of  freshly  mixed  agent,  the  newly  prepared  solution  was  left  undisturbed  for  some  minutes 
(5  to  10,  depending  upon  the  test  solution  concentration)  before  measurements  in  a  new  trial  were 
started. 

The  oscilloscope  was  configured  to  acquire  the  received  signal  and,  from  this  signal,  to 
generate  a  real  time  FFT  and  an  average  of  the  FFT  over  many  sweeps  (usually  100).  The  data  of 
interest  were  the  amplitude  peaks  of  the  fundamental  and  second  harmonic  frequencies  as  read  via 
the  cursors  from  the  averaged  FFT.  The  data  were  tabulated  for  serial  averages  over  the  duration 
of  the  run  to  record  the  change  in  the  amplitude  peaks  as  a  function  of  time. 

A  timer,  reset  with  the  initiation  of  each  average,  ensured  that  the  data  were  read  from  equal 
intervals  of  time;  the  intervals,  obviously,  were  longer  than  the  time  required  by  the  oscilloscope  to 
complete,  say,  a  100  sweep  average  plus  the  time  needed  to  read  the  data  from  the  display  (two 
minute  intervals  were  generally  used);  the  run  duration  was  usually  60  minutes. 

ni.2.2.  Contrast  agent  second  harmonic  characteristics 
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This  set  of  experiments  was  similar  in  method  to  that  described  in  the  preceding  section. 
Conceptually  the  major  difference  was  that,  whereas  in  the  previous  set  of  experiments  the 
characteristic  under  study  was  strictly  amplitude-versus-time,  this  set  of  experiments  sought  to 
define  the  contrast  agent  behavior  (represented  by  second  harmonic  and  fundamental  frequency 
amplitudes;  more  particularly,  the  ratio  of  these  amplitudes)  as  a  function  of  transmit  signal 
amplitude,  transmit  frequency,  contrast  agent  concentration  and  tone  burst  length.  The  results  of 
the  previous  experiments  guided  our  choice  of  the  length  of  the  data  runs  in  this  set  such  that  time 
was  not  a  significant  parameter. 

The  amplitude  ratio  characteristic  of  the  contrast  agent  may  be  perceived  as  occupying  a 
region  defined  by  a  hypercube  whose  axes  are  the  parameters  given  above.  Little  work  has  been 
published  regarding  the  characteristics  of  this  contrast  agent;  its  behavior,  in  this  context,  was  very 
poorly  understood.  Thus  it  became  essential  for  us  to  determine  the  practical  boundaries  of  the 
hypercube 

The  experimental  set  ups  consisted  of  both  the  two  transducer  configuration  described 
previously  (and  shown  in  figure  1)  and  the  single  transducer  set  up  of  figure  2,  which  will  be 
described  in  more  detail  here. 

One  transducer  was  installed  in  the  test  tank.  The  test  tank  could  be  configured  as  in  the 
previous  section;  collected  data  would  represent  only  the  reflection/scattering  of  second  harmonic 
energy  derived  from  fundamental  frequency  energy  impinging  axially  on  the  contrast  agent 
microbubbles  during  a  tone  burst  (the  single  transducer  configuration  implies  the  use  of  a  tone 
burst  drive).  Targets  could  also  be  installed  in  the  transducer  beam  to  provide  additional  data 
regarding  ranging  and  the  effect  of  interfacial  surfaces. 

The  active  (transmit/receive)  transducer  was  connected  to  the  common  RF  port  of  a  SPDT 
transmit/receive  switch.  One  throw  of  the  switch  connected  with  the  function  generator  output;  the 
other  throw  connected  to  the  oscilloscope  (either  with  or  without  intervening  preamplification). 

The  function  generator  gate  output  signal  provided  switching  control  synchronized  to  the  tone  burst 
output. 

Sets  of  experiments  were  categorized  by  the  parameter  under  investigation  (e.g.,  one  set  of 
experiments  consisted  of  runs  in  which  only  the  transmit  drive  voltage  was  varied,  all  other 
parameters  held  fixed).  Therefore,  data  were  collected  for  "drive  voltage"  experiments, 

"frequency"  experiments,  "concentration"  experiments  and  "tone  burst  length"  (i.e.,  the  number  of 
cycles  in  the  tone  burst)  experiments. 

Since  these  data  trials  were  much  shorter  than  those  in  the  lifetime  experiments,  a  set  of 
trials  could  be  completed  before  the  test  solution  required  renewal.  Thus,  for  a  given 
concentration,  a  trial  spanning  a  drive  voltage  interval  of  interest  could  be  repeated  over  a  span  of 
frequencies,  generating  a  two-dimensional  data  matrix.  Data  point  span  and  intervals  were 
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predetermined  (e.g.,  3  to  9  volts  of  drive  in  intervals  of  0.5  volt)  and  the  amplitude  information 
was  again  read  via  the  cursors  from  the  averaged  FFT  waveform. 

111.2.3.  Contrast  agent  shadowing  characteristics 

In  addition  to  the  reflection,  scattering  and  nonlinear  transformation  characteristics 
addressed  by  the  previous  experiments,  the  contrast  agent  exhibited  an  attenuation  or  absorption 
characteristic  which  is  commonly  known  as  "shadowing".  The  effect  of  shadowing  was  to  reduce 
the  returned  signal  by  up  to  100%,  depending  on  the  environment.  The  effect  was  noticeable  when 
targets  were  used  and  when  they  were  not  used  (in  effect,  the  contrast  agent  shadowed  itself).  The 
mechanism  causing  this  effect  is  well  known  and  related  to  the  bubbly  nature  of  the  agent;  bubble 
are  highly  efficient  reflectors  of  acoustic  energy,  and,  when  present  in  high  enough  concentration, 
will  obscure  all  energies  distal  to  the  incident  beam.  Incidental  results  from  other  experiments 
convinced  us  that  it  would  be  worthwhile  to  attempt  characterization  of  this  behavior. 

The  experimental  set  up  was  predominately  that  illustrated  in  Figure  5,  with  targets  in  the 
test  tank.  The  characterization  was,  again,  amplitude  versus  time.  The  differences  from  the 
lifetime  experiments  were:  1)  amplitudes  of  the  averaged  real  time  received  signal,  not  of  the 
averaged  FFT  were  measured;  2)  rather  than  a  decay  in  amplitude  with  time,  we  were  looking  for  a 
decay  in  shadowing.  The  aim  was  to  determine  the  concentration  of  the  contrast  agent  based  on 
amplitude  readings  of  received  signals  reflected  from  axially  aligned  targets. 

An  experimental  trial  consisted  of  serial  amplitude  readings,  from  at  least  two  targets, 
starting  10  seconds  after  the  introduction  of  the  contrast  agent  into  the  tank  and  continuing  until  a 
steady  state  condition  became  established.  The  readings  were  performed  as  rapidly  as  possible, 
utilizing  the  oscilloscope's  cursors,  and  noting  the  time  each  reading  occurred. 

111.2.4.  Ex  vivo  experiments 

Figure  7  illustrates  an  experimental  set  up  which  utilizes  a  biological  preparation,  in  this 
case  an  excised  rabbit  kidney.  The  kidney  was  immersed  (after  being  purged  with  normal  saline 
injected  into  the  renal  artery)  into  a  solution  of  normal  saline  (0.9%  NaCl).  Small  bore,  flexible 
tubing  was  then  cannulated  to  both  the  renal  artery  and  vein.  The  venous  line  lead  to  another 
plastic  container  which  served  as  a  venous  reservoir  to  increase  the  circulating  volume  to  a  more 
physiologic  amount  (so  that  the  contrast  agent  injections  would  experience  a  more  realistic  dilution 
while  circulating)  and  to  accommodate  volume  changes  caused  by  serial  agent  injections.  Another 
flexible  line  connected  the  reservoir  to  the  inlet  port  of  the  circulating  pump.  The  bellows  pump 
head  provided  a  rough  approximation  of  a  pulsatile  outflow  through  the  flexible  plastic  tubing 
cannulated  to  the  renal  artery.  On  the  tubing,  between  the  pump  outlet  and  the  renal  artery,  was 
installed  an  injection  port,  to  which  was  attached  a  three-way  stopcock.  On  the  remaining  two 
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ports  of  the  stopcock  were  attached  syringes  which  contained  contrast  medium  and  a  saline  flush 
solution. 

The  contrast  agent  was  prepared  by  diluting  a  weighed  amount  of  the  granules  (generally 
between  300  and  900  mg)  with  a  small  volume  of  normal  saline,  in  a  glass  vial  to  achieve  the 
desired  concentration  (generally  150  to  450  mg/ml).  The  solution  was  then  withdrawn  into  a 
syringe  and  attached  to  the  stopcock  in  preparation  for  injection.  The  contrast  medium  was  injected 
as  a  bolus,  followed  by  a  small  bolus  injection  of  flush  solution. 

The  transducer,  supported  on  the  arm  of  a  micromanipulator  assembly  (not  shown)  had 
been  positioned  over  the  target  region  of  the  kidney  before  contrast  injection;  the  electronics  set  up 
was  similar  to  that  of  Figure  5.  Data  acquisition  began  with  the  completion  of  the  flush  injection 
and  continued  until  the  signal  amplitude  decreased,  because  of  dilution,  to  the  noise  floor,  or  until  a 
steady-state  response  was  noted.  Again,  amplitude  peaks  were  read  via  the  display  cursors  from 
either  the  averaged  FFT  (for  second  harmonic  studies)  or  from  the  averaged  real  time  received 
signal  (for  shadowing  studies). 

IV.  EXPERIMENTAL  RESULTS 
IV.  1.  Lifetime  experiments 

The  lifetime  data  were  to  provide  information  we  believed  would  illuminate  such  properties 
of  the  contrast  agent  as  related  to  subsequent  experimental  protocols.  We  wanted  to  know  how  the 
agent  changed  with  time,  how  repeatable  the  data  from  a  given  concentration  would  be,  how 
handling  might  affect  the  agent's  properties  and  generally  how  the  diluted  agent  behaved  in  contact 
with  the  experimental  environment. 

Figure  8  is  a  graph  of  a  data  set  typical  of  this  group  of  experiments.  The  amplitude  (in 
mV)  vs.  time  (in  minutes)  behavior  of  the  fundamental  frequency  peak,  the  second  harmonic 
frequency  peak  and  the  ratio  (second/fundamental)  are  plotted.  Characteristic  of  this  group  of 
experiments  was  the  relatively  constant  value  of  the  amplitude  peak  ratios  during  a  trial.  For  the 
trial  depicted  in  Figure  8  the  average  ratio  was  14.3%  with  a  standard  deviation  of  2.35%.  The 
ratio  exhibited  a  slight  positive  correlation  with  time  (average  ratio  13.92%  during  the  first  20 
minutes  of  the  run;  14.6%  during  the  last  20  minutes)  while  the  standard  deviation  of  the  ratio 
presented  with  a  slightly  negative  correlation  (2.73%  in  the  first  20  minutes;  2.02%  for  the  last 
20).  The  data  bore  out  our  subjective  impression  of  a  high  degree  of  kinetic  activity  in  a  newly 
mixed  solution,  which  decreases  with  time.  Stirring  the  test  solution  during  the  course  of  a  run 
seemed  to  produce  a  subtle  effect  on  the  ratio  (inferring,  of  course,  an  effect  on  the  fundamental 
and  second  harmonic  amplitude  peaks),  as  shown  in  Figure  9.  Concentration  had  little  effect  on 
the  characteristic  profile  of  the  time-dependent  behavior. 
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During  this  series  of  experiments,  we  also  found  that  the  plastic  film  septum  (a  technique 
carried  over  from  earlier  work)  used  on  the  transducer  holding  fixtures  was  unnecessary  in  the 
configuration  employed  in  the  in  vitro  experiments  of  this  project.  Removal  of  the  film  also  had 
the  salutary  effect  of  increasing  the  sensitivities  of  the  acoustic  elements  of  the  system. 

IV.2.  Second  harmonic  experiments 

We  expected  this  group  of  experiments  to  provide  the  information  needed  to  determine 
which  portion  (or  portions)  of  the  hypercube  were  most  suitable  for  this  project.  Much  of  the 
effort  in  these  experiments  went  toward  exploring  the  extremes  of  the  parameters  of  interest,  i.e., 
the  comers  and  edges  of  the  hypercube.  Correspondingly,  drive  levels  over  a  range  of  200  mV 
(for  CW)  to  30  V  (for  tone  burst),  concentrations  of  0.25  mg/ml  to  5  mg/ml,  burst  lengths  of  2  to 
12  cycles  and  frequencies  of  1.5  to  3  MHz  were  utilized. 

Figure  10  is  a  graph  of  representative  data  from  this  family  of  experiments.  The  graph 
depicts  the  behavior  of  the  second  harmonic/fundamental  peak  amplitude  ratio  as  a  function  of 
transmit  drive  voltage  with  the  tone  burst  length  as  a  parameter.  The  contrast  agent  concentration 
was  1  mg/ml  and  the  drive  frequency  was  2.5  MHz.  Each  data  point  was  obtained  from  the 
averaged  FFT  trace  on  the  oscilloscope,  an  averaged  trace  consisting  of  250  sweeps.  The 
experimental  set  up  was  that  of  Figure  4. 

The  amplitudes  of  the  fundamental  and  the  second  harmonic  generally  shifted  upward,  at  a 
given  drive  level,  as  the  tone  burst  length  increased  (see  Figures  10a,  b,  c,  d)  and  the  ratio  shifted 
down.  Although  gross  trends  in  the  curves  could  be  related  (e.g.,  both  the  fundamental  and 
second  harmonic  amplitudes  roughly  increased  with  drive  voltage,  at  a  given  tone  burst  length) 
from  one  run  to  the  other,  it  is  clear  that  details  of  the  curves  are  significantly  variable.  Only  one  of 
the  curves  is  monotonic. 

Not  unexpectedly,  we  found  that  CW  experiments  required  a  lower  transmit  drive 
amplitude  than  tone  burst  experiments  and  tolerated  a  lower  contrast  concentration  while  still 
providing  useful  data  (i.e.,  the  frequency  peaks  were  easily  distinguished  from  the  FFT  noise 
floor). 

IV. 3.  Shadowing  experiments 

Figure  1 1  is  a  plot  of  a  screen  dump  from  the  LeCroy  oscilloscope  (to  an  HP  7470  plotter). 
The  experimental  set  up  was  that  depicted  in  figure  1  with  the  addition  of  the  targets  shown  in 
Figure  5.  The  top  trace  in  Figure  1 1  is  the  receive  signal;  the  echoes  from  each  target  are  easily 
distinguished.  Though  the  second  echo  has  traversed  a  greater  path  length,  the  target  which  caused 
that  echo  presented  a  better  reflecting  surface  than  did  the  nearer  target;  thus  the  greater  amplitude 
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in  the  echo  from  the  farther  target.  The  origin  of  the  third  echo  was  probably  due  to  reflection  of  a 
high  amplitude  signal  from  the  walls  of  the  test  chamber. 

Averaging  the  receive  signal  over  50  sweeps  produced  the  bottom  trace;  the  baseline 
acoustic  and  electronic  noise  is  significantly  reduced.  Averaging  also  minimized  the  effects  of  the 
slight  amplitude  jitter  present  on  the  receive  signal.  The  step  displacements  in  the  trace  2  baseline 
are  an  artifact  of  the  oscilloscope. 

Asymmetry  of  the  echo  waveforms  about  the  baseline  precluded  reading  simply  the  peak 
values  of  each  echo.  Maximum  and  minimum  amplitudes  (which  provided  the  peak-to-peak 
amplitude)  for  each  target  echo  were  read  via  the  cursors. 

Figure  12  is  a  graph  of  representative  data  from  this  series  of  experiments  (note  the  change 
in  time  interval  for  the  last  two  data  points).  The  salient  feature  of  this  data  plot  is  the  difference  in 
attenuation  (due  to  the  difference  in  path  length)  of  the  return  echoes  from  targets  1  and  2.  The 
theoretical  development  (see  appendix)  for  utilizing  the  shadowing  characteristic  required  that 
attenuation  be  a  function  of  path  length;  these  experiments  verified  the  theoretical  treatment. 

IV.4.  Ex  vivo  experiments 

The  bench  top  in  vitro  experiments  provided  practical  experience  with  the  contrast  agent. 

To  better  approximate  the  intended  environment  for  the  application,  however,  experiments  utilizing 
a  biological  preparation  were  important  to  assess  the  validity  of  applying  in  vitro  results  to  our 
conclusions  regarding  the  feasibility  of  the  proposed  technique. 

Attempts  to  elicit  a  second  harmonic  response  were  equivocal  in  this  set  of  experiments. 
We  were  unable  to  record  any  reliably  quantifiable  second  harmonic  data.  Given  past  experience 
with  the  agent  in  Doppler  experiments,  we  propose  that  this  was  due  to  the  insensitivity  of  our 
transducer  system  with  regard  to  the  second  harmonic  in  a  highly  scattering,  biologic  environment. 
(We  did  attempt  to  obtain  a  more  sensitive  probe;  time  limitations,  however,  made  it  impossible  to 
integrate  it  and  test  it  in  the  system  before  the  end  of  the  Phase  I  effort.) 

The  shadowing  phenomenon  was  relatively  easy  to  explore  with  the  rabbit  kidney, 
however.  Figure  13  is  a  graph  of  the  two  trials  in  which  we  specifically  investigated  the 
attenuation  characteristic  of  the  contrast  agent.  The  trials  were  separated  by  about  15  minutes,  and 
demonstrate  encouraging  repeatability.  The  data  from  trial  #2  presents  a  generally  flatter  profile 
with  slightly  less  attenuation.  In  both  trials,  steady  state  conditions  occurred  approximately  2 
minutes  after  injection.  The  dips  in  the  curves,  we  believe,  result  from  recirculation  of  the  contrast 
agent.  The  agent  possesses  a  long  half  life  and  the  background  level  contributed  from  trial  #1 
explains  the  previously  noted  characteristic  of  the  trial  #2  curve. 

Figure  14  is  a  graph  of  a  trial  which  is  representative  of  a  series  performed  during  a 
subsequent  experiment  (in  this  case,  the  ordinate  is  attenuation  rather  than  amplitude).  Again, 
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steady  state  occurs  at  approximately  2  minutes  post-injection,  the  effect  of  residual  agent  from 
preceding  injections  is  evident  during  serial  trials,  recirculation  is  variably  evident  from  run  to  run 
and  the  general  profile  of  the  attenuation  vs.  time  behavior  of  the  contrast  agent  is  very  similar 
among  the  trials. 

Several  of  the  trials  performed  during  this  group  of  experiments  explored  the  effects  of 
concentration.  As  anticipated,  attenuation  displayed  a  positive  correlation  to  concentration  of  the 
injectate.  Figure  15  depicts  the  results  of  an  in  vitro  experiment  to  also  determine  the  effect  of 
contrast  agent  concentration  on  shadowing.  The  curves  resulting  from  the  concentration  runs  are 
easily  distinguished  from  each  other. 

V.  CONCLUSIONS 

All  objectives  from  the  proposal  have  been  met.  Additional  comments  regarding  the 
Soldier  Individual  Computer  program  and  results  from  experiments  conducted  during  this  project 
suggest  the  form  of  future  work. 

Foremost  is  the  question  of  appropriateness  of  this  technique  relative  to  the  Soldier 
Individual  Computer  project.  The  methods  under  consideration  here  (second  harmonic  and 
shadowing)  require  injection  of  a  contrast  agent.  Although  this  precludes  the  possibility  of 
continuous  monitoring  of  soldier  health,  conceivably  a  baseline  data  could  be  obtained  under 
controlled  conditions  and  stored  and  used  to  compare  with  data  from  an  injury  or  trauma  situation, 
at  which  time  another  injection  would  be  made. 

The  second  harmonic  technique  has  shown  significant  promise  of  practical  application  in 
this  and  other  works,  and  with  improvements  in  understanding,  technique  and  technology  could 
develop  into  a  viable  clinical  approach  to  quantifiable  measurement  of  blood  flow.  Thus,  future 
efforts  will  continue  to  be  directed  toward  rigorous  testing  of  the  agent,  especially  in  in  vivo 
conditions.  Rapid  advances  in  computing  power  will  obviously  aid  in  development  of  a  device  (in 
this  case,  a  computer  chip  for  the  Soldier  Individual  Computer). 

The  well  known  phenomenon  of  shadowing  demonstrates  the  immediately  attractive 
attributes  of  predictability  and  docility.  However,  until  now,  this  has  been  viewed  as  an 
annoyance  to  be  avoided  in  ultrasound  diagnostics,  never  studied  for  practical  use.  The 
preliminary  work  performed  here  shows  much  promise  for  perfusion  quantification,  and  further 
work  will  elaborate  upon  these  experiments. 
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VII.  APPENDIX  1 :  Theoretical  Development  of  Shadowing  Phenomenon 

In  addition,  we  have  serendipitously  discovered  an  interesting  phenomenon  related  to  the 
prominent  shadowing  effect  of  the  contrast  agent.  This  is  a  repeatable  effect,  and  seems  to  have 
stable  behavior.  Based  on  preliminary  observations  from  this  effect,  we  have  developed  a  new 
theory  to  utilize  this  phenomenon  for  perfusion  quantification.  This  will  briefly  be  described  here; 
a  full  analysis  will  of  course  be  presented  in  the  final  project  report. 

Assuming  that  the  backscattered  signal  from  a  target  located  at  a  distance  x  from  a  pulsed 
system  can  be  described  by: 

r(x)  =  TpG(x)e  2ax  (VII  1) 

where  r  is  the  received  signal,  T  is  the  transmitted  signal,  p  is  the  target  reflection  coefficient,  a  is 
a  constant  attenuation  coefficient  and  G  is  the  transducer  pattern  envelope  function.  If  we  assume 
that  the  attenuation  coefficient  is  varying  across  the  medium,  equation  (VII.  1)  takes  the  following 
form: 


r(x)  =  TG(x)pe  j 


2ct(x)dx 


(VII.  2) 
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If  we  assume  that  the  attenuation  coefficient  in  a  small  region  can  be  considered  constant,  we  can 
write  the  following  two  equations  for  the  received  echoes  from  the  targets.  The  received  echo  from 
the  first  target  would  be: 


rx  i 


r(xi)  =  TG(xi)p!e'J0 
and  from  the  second  target: 


2a(x)dx 


(VII.  3) 


r(x2)  =  TG(x2)p2e  J0 


2a(x)dx 


=  TG(x2)p2e  J0 


2a(x)dx  -  2a(x)dx 


=  TG(x2)p2e 


-2axAx  -  2a(x)dx 


(VII.4) 

where  ax  is  the  "constant"  attenuation  coefficient  between  the  two  targets  and  Ax  is  the  distance 
between  the  targets. 

Now,  assuming  that  the  marginal  attenuation  coefficient  caused  by  a  perfused  contrast 
agent  is  directly  proportional  to  the  contrast  agent  density  (concentration).  Thus,  the  marginal 
attenuation  coefficient  due  to  the  contrast  agent  will  be  an  =  yN.  So,  in  general,  we  add  the  term 
a(x)  to  the  term  yN(x)  to  estimate  the  total  attenuation  coefficient. 

Equations  (VII.3)  and  (VII.4)  describe  the  received  signal  from  two  closely  spaced  targets 
before  contrast  agent  injection.  After  contrast  agent  injection  and  with  the  added  assumption  that 
the  contrast  agent  density  between  the  two  targets  can  be  assumed  constant,  equations  (VII.3)  and 
(VII.4)  become: 


rn(xi)  =  TG(xi)pie 


2(yN(x)  +  a(x))dx 


(VII.5) 


rn(x2)  =  TG(x2)p2e 


-2(ax+yNx)Ax- 


2(tN(x)  +  a(x))dx 


Jo 


(VII.  6) 

Denoting  the  ratio  between  equation  (VII.4)  and  equation  (VII.3)  as  R  and  between  equations 
(VII.6)  and  (VII.5)  as  Rn  we  obtain: 
p2G(x2)e-2a>‘Ax 


R  = 


PiG(xi) 


(VII.7) 


_  p2G(x2)e-2(a*+'yN*)Ax 
piG(xi) 


(VII.  8) 


From  equations  (VII.7)  and  (VII.8)  we  can  extract  Nx: 


Nx 


(VII.9) 


2yAx 

Equation  (VII.9)  represents  the  solution  of  density  of  contrast  agent,  which  is  presumed  is  directly 
proportional  to  amount  of  tissue  perfused  (between  the  two  targets).  Further  refinement  of  this 
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theory  to  in  vivo  situations,  as  well  as  further  experimentation  to  verify  the  assumptions,  is  in 
progress. 

Vm.  APPENDIX  2:  Theoretical  Development  of  Second  Harmonic  Quantification 

It  has  been  shown  in  previous  analyses  that  will  not  be  repeated  here  (Schrope  et  al,  1992; 
Schrope  and  Newhouse,  1993)  that  ultrasound  by  itself  is  unsuitable  for  detection  of 
microcirculatory  blood  flow,  especially  that  located  in  deep  tissues.  When  a  specific  ultrasound 
contrast  agent  (SHU-508  or  Levovist®,  Schering  AG  and  Berlex  Laboratories)  is  added  to  the 
system,  however,  a  unique  "signature"  is  applied  to  the  slow,  small  signal  blood  flow  that  allows  it 
to  be  unambiguously  detected  over  other  noise  in  the  system,  including  signal  from  surrounding 
tissue  movements.  The  contrast  agent,  consisting  of  a  suspension  of  microbubbles  on  the  order  of 
two  to  three  microns  in  diameter,  is  administered  in  a  small  volume  (usually  less  than  5  ml) 
intravenous  bolus.  In  the  present  application,  the  three-step  management  of  critical  trauma  victims, 
that  is,  Airway  — >  Breathing  — >  Circulation,  dictates  that  an  intravenous  line  be  placed  for 
administration  of  various  medications  and  fluids  as  needed.  This  IV  site  can  also  be  used  to 
administer  the  agent  necessary  for  the  measurement.  Further  as  yet  unpublished  work  is  presented 
here  to  describe  the  efforts  at  quantification  of  the  data  as  a  measure  of  blood  perfusion. 
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Figure  VIII 1:  Signals  and  attenuation  coefficients  used  in  the  following  analysis  and  should  be  used  as  a  reference 
block. 


Figure  VIII1  shows  signals  and  attenuation  coefficients  used  in  the  following  analysis  and 
should  be  used  as  a  reference  block  diagram.  We  assume  the  first  and  second  harmonic  signal 
voltages  (in  the  analysis,  we  use  voltages  and  pressures  interchangeably)  of  the  backscattered 
signals  near  the  bubble  site  can  be  represented  by  the  following  expressions: 

vi  =  wo  (VIII.  1) 
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(VIII.2) 


v2  =  dam  +  a2v%) 

According  to  eq.  (VIII.  1),  the  first  harmonic  backscattered  pressure  signal  v\  near  the  bubble  site 
is  directly  proportional  to  the  transmitted  signal  strength  near  the  bubble  site  multiplied  by  the 
scattering  cross  section  of  the  bubble,  Cp  Likewise,  the  second  harmonic  signal  is  proportional  to 
the  bubble  scattering  cross  section  multiplied  by  a  second  order  polynomial  expression. 

The  transmitted  signal  Vo  is  attenuated  by  frequency  dependent  attenuation  of  the  tissue 
which  is  characterized  by  the  coefficient  Kt(r)  (for  clarity,  we  omit  the  dependence  on  r  in  the 
following  analysis  and  assume  a  set  range  delay)  and  by  transducer-skin  and  organ  wall  interface 
losses  which  are  represented  by  the  coefficient  Kx.  The  signal  strength  near  the  bubble  site  is  also 
determined  by  the  transducer  beam  geometry  Kw.  vo  can  thus  be  written  as: 

vo  =  V0KxlKnKwl  (VIII.3) 

The  1  added  to  the  subscripts  denotes  that  the  attenuation  and  beam  geometry  coefficients  apply  to 
the  first  harmonic  frequency.  The  signal  vi  (the  first  harmonic  backscattered  pressure  generated  by 
a  single  bubble)  is  according  to  eq.  (VIII.  1): 

V!  =  oV  qK  X\K  t\Kw\  (VIII.4) 

The  average  backscattered  voltage  from  an  ensemble  of  N  randomly  distributed  bubbles  at  the 
receiving  transducer  is: 

Vi  =  VoWaK^K^  (VIII.5) 

The  backscattered  voltage  is  directly  proportional  to  the  square  root  of  the  number  of  bubbles  in  the 
range  cell  (as  the  received  power  is  directly  proportional  to  N).  Reflecting  back  to  the  transducer, 
the  backscattered  signal  is  attenuated  again  by  the  Kxi,  Kti  and  Kwi  coefficients  described  above. 

For  the  second  harmonic  signal  for  a  single  bubble  (see  eq.  (VIII.2)) 

v2  =  dam  +  a2v%) 

v2  =  daiV0KxlKtlKwl  +  a2VlK2xlK}lKll)  (VIII.  6) 

and  the  signal  received  by  the  transducer  (second  harmonic)  is: 

L2  =  v2Kx2Kt2Kw2  (VIII.7) 

The  2  added  to  the  subscripts  denotes  that  the  attenuation  and  beam  geometry  coefficients  apply  to 
the  second  harmonic  frequency.  Note  that  the  attenuation  and  transducer  beam  coefficients  for  the 
second  harmonic  signal  are  different  than  for  the  first  harmonic  signal.  The  received  second 
harmonic  signal  is  thus  determined  by  both  the  first  harmonic  coefficient  (transmitted)  and  the 
second  harmonic  coefficients  (reflected) 

V2  =  daiV()KxlKnKwl  +  a^oK^K^K^^K  x2K  t2K  w2  (VIII.8) 

The  received  signal  from  N  bubbles  is  thus: 

V2  =  V()fN  oK X\K t\K W\K x2K t2K w^a\  +  a2V qK X\K t\K w\]  (VIII.9) 
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The  result  of  eq.  (VIII.9)  represents  the  general  expression  for  the  second  harmonic  received  signal 
from  an  ensemble  of  randomly  distributed  scatterers.  We  will  use  this  fundamental  result  in  the 
next  section  to  derive  the  attenuation  expression  and  then  solve  for  the  number  N  of  the  scatterers 
in  the  range  cell. 


Attenuation  Estimation 


The  attenuation  can  be  derived  from  the  ratio  of  the  received  second  harmonic  signals  for 


two  different  transmitted  signals,  say  Vqi,  Vq2-  Denote  the  received  second  harmonic  signals 


respectively  as  V21,  V22-  We  get  (from  eq.  (VIII.9)): 

V21  _  Voitfll  +  a2KX\Kt\Kw\Vo\\ 

V22  Voia\  +  a2^ X\K  t\Kw\V  Q2\ 


(VIII.  10) 


Denote: 


P  =  KxlKnKwl 

n  -  ^21.  n  _  VQI 

T  7  ) 

y  22  VQ2 


(attenuation  to  bubble  -  first  harmonic) 


(v1n.1i) 


where  V21,  V22  are  the  second  harmonic  received  signals  generated  by  the  contrast  agent  due  to 
V01,  V02  respectively.  We  get  from  eq.  (VIII.10)  and  eq.  (VIII.  1 1): 

Ro 


1+^pVoi 


and 


<22(^0  } 

P  _  ^t\^< 

Voi-^V02 


(VIII.  12) 


(VIII.  13) 


If  a2/ai  is  given,  the  number  of  bubbles  in  the  range  cell  can  be  extracted  using  eq.  (VIII.5) 
Vn  =  oVoiWp2 

where  Vn  is  the  first  harmonic  received  signal  due  to  the  transmitted  signal  Vqi. 


N 


V 


11 


(VIII.  14) 


The  challenge  is  to  obtain  Vn,  i.e.,  to  subtract  the  tissue  echoes  from  the  received  first  harmonic 
signal  that  contains  both  the  bubbles  as  well  as  tissue  echoes. 


Calibration 

In  the  previous  section,  we  assumed  that  the  ratio  ai/a2  was  known.  In  some  cases  this 
value  may  not  be  provided  by  the  contrast  agent  vendor  and  it  must  be  evaluated.  This  can  be  done 
in  the  following  manner:  in  an  experiment  where  N,  a,  Vo  and  V\  are  known,  the  attenuation  term 
(3  can  be  extracted  using  eq.  (VIII.5) 
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<VI11-15) 

ai/a2  can  now  be  calculated  using  eq.  (VIII.  15)  using  two  different  transmitted  signals.  In  order  to 
improve  the  accuracy  of  the  quantity  (3,  multiple  experiments  with  a  variety  of  transmitted  signal 
voltages  Vo  as  well  as  bubble  concentration  N  should  be  carried  out  and  averaged.  Likewise,  to 
improve  the  estimate  accuracy  of  the  ratio  ai/a2  several  experiments  with  different  transmitted 
signal  voltages  should  be  carried  out  and  averaged.  This  should  also  provide  a  test  to  the  second 
harmonic  signal  generation  model  that  assumes  that  two  terms  in  the  polynomial  expansion  of  eq. 
(VIII.2)  are  sufficient.  If  the  ratio  ai/a2  shows  a  significant  trend  (instead  of  just  random 
experimental  variations),  the  model  should  probably  be  refined  as  will  be  discussed  later. 

A  practical  way  to  calibrate  the  system  can  be  achieved  by  mixing  a  known  quantity  of  the 
contrast  agent  with  a  known  quantity  of  water  to  create  a  "standard"  to  which  all  other 
measurements  will  be  related.  In  this  case  we  first  estimate  (3  using  eq.  (VIII.  15).  Note  that  the 
quantity  defined  by  this  "standard"  is  of  No2,  and  not  only  of  concentration  N.  This  standard  does 
not  have  to  be  a  close  approximation  of  the  actual  quantity  of  No2.  It  can  be  set  arbitrarily  to  1,  for 
example.  This,  naturally,  will  cause  the  ratio  ai/a2  to  be  tied  to  this  arbitrary  "standard".  While 
not  affecting  the  relative  quality  of  concentration  estimates,  it  may  cause  confusion  if  a  standardized 
"standard"  is  not  defined. 

Microvctscular  Blood  Flow  Estimation 

In  this  section  we  describe  the  underlying  fundamental  theory  of  the  proposed  innovation. 
While  there  have  been  quite  a  few  attempts  in  the  past  to  use  ultrasound  for  microcirculatory  blood 
flow  estimation,  none  of  the  proposed  techniques  provides  adequate  solution  to  the  problem  (see 
section  c). 

The  main  problem  with  the  current  methods  is  that  the  transmitted  and  received  signals  are 
attenuated  by  unknown  amount,  by  highly  varying  tissue  media,  blood  vessel  walls  and 
transducer-skin  interface.  It  is  thus  impossible  to  estimate  accurately  the  reflection  from  the  tissue, 
even  though  contrast  agents  for  image  contrast  enhancement  are  available. 

Our  breakthrough  approach,  as  was  explained  earlier,  is  based  on  the  use  of  a  unique 
contrast  agent,  SHU-508  (Schering  AG  and  Berlex  Laboratories).  However,  unlike  other 
approaches  using  contrast  agents,  such  as  area  under  the  curve  or  half-lifetime  (Monaghan  et  al, 
1988;  Ong  et  al,  1984;  Ten  Cate  et  al,  1984),  our  method  is  based  on  the  nonlinear  characteristics 
of  this  particular  contrast  agent.  These  agents  that  typically  look  like  spheres  produce  in  many 
cases  relatively  strong  signals  at  twice  the  frequency  of  the  transmitted  signal  (Schrope  et  al,  1992; 
Schrope  and  Newhouse,  1993;  Bums  et  al,  1993).  We  refer  to  these  signals  as  second  harmonic 
signals.  In  the  proposed  method,  an  ultrasonic  signal  with  a  given  frequency  is  transmitted 
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towards  the  region  of  interest  (which  is  the  region  for  which  we  attempt  to  estimate  the 
microcirculatory  blood  flow).  The  contrast  agent,  due  to  its  nonlinear  nature,  generates  a  second 
harmonic  signal  which  is  received  by  the  receiving  transducer.  The  first  harmonic  signal  is  also 
received  by  a  receiving  transducer  (possibly  the  same  transducer). 

The  basic  idea  behind  the  proposed  technique  is  that  the  characteristics  of  the  medium,  i.e., 
its  attenuation  characteristics,  can  be  extracted  from  the  ratio  of  the  received  second  harmonic 
signals  from  two  different  transmitted  signal  amplitudes.  For  example,  we  transmit  first  a  signal  at 
10  V  and  record  the  received  second  harmonic  signal  amplitude  from  a  region  of  interest,  say,  10 
mV.  We  repeat  this  experiment  with  a  different  transmitted  signal,  say  20  V,  and  record  the 
received  second  harmonic  signal  amplitude  for  this  case,  say  30  mV.  Note  that  due  to  the 
nonlinear  behavior  of  the  contrast  agent  the  second  harmonic  tripled  for  only  a  two-fold  increase  in 
the  transmitted  signal  strength.  As  the  following  analysis  will  clearly  show,  tissue  attenuation 
coefficients  can  be  extracted  from  the  ratio  of  the  two  second  harmonic  signals  with  respect  to  the 
ratio  of  the  transmitted  signals  (in  this  case  the  output  ratio  is  3  =  30/10,  while  the  input  ratio  is  2  = 
20/10).  The  main  reason  that  we  can  estimate  the  attenuation  is  due  to  the  fact  that  the  ratio 
between  the  two  second  harmonic  signals  depends  solely  on  the  strength  of  the  transmitted  signal 
at  the  contrast  agent  site.  The  stronger  the  transmitted  signals  at  the  contrast  agent  site,  the  larger 
the  second  harmonic  signal  ratio  would  be.  Thus,  we  can  infer  from  this  ratio  the  signal  strength  at 
the  contrast  agent  site  and  consequently  the  medium  attenuation.  In  practice,  the  two  different 
transmitted  signals  would  be  transmitted  sequentially  in  a  rapid  rate,  say  every  0.5  to  1  msec.  This 
way  the  spheres  that  compose  the  contrast  agent  do  not  change  significantly  their  random  relative 
positions  and  the  target  can  be  considered  to  be  unchanged. 

Once  the  attenuation  has  been  quantified,  we  would  "isolate"  the  echo  from  the  contrast 
agent  and  estimate  very  accurately  the  concentration  of  the  contrast  agent.  As  was  explained 
earlier,  the  contrast  agent  is  injected  in  a  controlled  fashion.  Thus,  the  concentration  of  the  injected 
contrast  agent  is  known.  As  a  consequence  an  accurate  estimate  of  blood  perfusion  would  be 
possible. 
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Figure  1:  Block  diagram  of  system  design  number  one,  based  on  existing  off-the-shelf 
components. 
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range-delay 


Figure  2:  Block  diagram  of  system  design  number  two,  where  much  of  the  processing  has  been 
transformed  to  DSP. 


Figure  3:  Block  diagram  of  system  design  number  three,  with  Doppler  capability. 
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Figure  4:  In  vitro  experimental  apparatus  -  two  unfocused  transducers,  tone  burst  transmitted 
signal. 
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Figure  5:  In  vitro  experimental  apparatus  -  single  transducer,  tone  burst  transmitted  signal. 
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Transducer  cable 


Figure  6:  Detail  of  transducer  fitting  assembly. 
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Figure  7:  Ex  vivo  experimental  preparation,  single  transducer  set  up 


DAMD17-95-C-5043;  Final  Report 


28 


Voltage/Ratio 


Figure  8:  Received  signal  amplitude  versus  time  for  lifetime  experiments,  unstirred  condition. 
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P^  -  echo  from  needle  target 

^2  =  echo  from  aluminum  strip  target 
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Figure  11:  Oscilloscope  screen  dump  of  target  echoes,  shadowing  experiment. 


Figure  12:  In  vitro  shadowing  experiment  -  concentration  of  contrast  agent  =  0.25  mg/ml,  transmit 
frequency  =  3.59  MHz,  drive  voltage  =  4.8  volts  (peak-to-peak),  burst  length  =  2  cycles;  data 
averaged  50  times. 


RETURNED  SIGNAL  FROM  A  RABBIT  KIDNEY 

350  mg /ml  1m!  injected  2.25  MHz 


Figure  13:  Ex  vivo  shadowing  experiment  -  time-dependent  characteristic  of  raw  received  signal. 
The  two  lines  depicted  indicate  two  different  trials. 
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Figure  14:  Shadowing  (or  attenuation)  as  a  function  of  time.  Contrast  agent  concentration  =  300 
mg/ml,  injected  volume  =  1  ml. 


LN(R/Rn)/2D  as  a  function  of  time  for 
three  levels  of  agent  concentration 
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Figure  15:  In  vitro  shadowing  experiment 
LN(R/Rn)/2D  parameter. 


see  theoretical  development  for  explanation  of 
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